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Abstract 



We present here new transmission spectra of the hot Jupiter HD-189733b 
using the SpeX instrument on the NASA Infrared Telescope Facility. We ob- 
tained two nights of observations where we recorded the primary transit of the 
planet in the J-, H- and K-bands simultaneously, covering a total spectral range 
from 0.94 to 2.4/zm. We used Fourier analysis and other de-trending techniques 
validated previously on other datasets to clean the data. We tested the statis- 
tical significance of our results by calculating the auto-correlation function, and 
we found that, after the detrending, white noise dominates at most frequencies. 
Additionally, we repeated our analysis on the out-of-transit data only, showing 
that the residual telluric contamination is well within the error bars. While these 
techniques are very efficient when multiple nights of observations are combined 
together, our results prove that even one good night of observations is enough 
to provide statistically meaningful data, which might appear counterintuitive 
given the daunting accuracy to be achieved. Our observed spectra are consis- 
tent with space-based data recorded in the same wavelength interval by multiple 
intruments, indicating that ground-based facilities are becoming a viable and 
complementary option to spaceborne observatories. The best fit to the features 
in our data was obtained with water vapour. Our error bars are not small enough 
to address the presence of additional molecules. 

Subject headings: techniques: spectroscopic, methods: data analysis, planets and 
satellites: atmospheres, planets and satellites: individual(HD- 189733b) 
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1. Introduction 

With over 800 exoplanets discovered and more to come in the near future, exoplanet 
research is today anything but dull. The synergy between space and ground work is 
delivering outstanding results and the menagerie of exoplanets is getting more diverse. 
For instance the Kepler mission (Borucki et al. 2010, 2011) has delivered, among other 
impressive discoveries, the smallest exoplanet known to date, which is Mars-sized (Muirhead 
et al. 2012), and the high- resolution spectrograph HARPS at La Silla Observatory 
(European Southern Observatory) recently detected a new Earth-sized exoplanet around a 
Centauri B (Dumusque et al. 2012). 

Beyond the key parameters of mass, radius and orbital inclination, the next critical 
step is to determine the chemical composition of these exotic worlds. Detecting atmospheric 
features, which have a contrast of about 10 -4 compared to the host star's radiation, is quite 
a challenge. However, for transiting exoplanets this has proved to be feasible from space 
and the ground. With the high stability of the Spitzer Space Telescope and the Hubble 
Space Telescope, the spectra of bright close-in massive planets have been obtained and 
ionic/atomic/molecular species like ionized forms of hydrogen, silicon, carbon, magnesium, 
as well as alkali metals, water vapor, methane, carbon monoxide and dioxide have been 
detected in a handful of these planets. With the demise of "cold Spitzer" and NICMOS, 
though, the wavelength range covered from space has narrowed. On the positive side the 
new Wide Field Camera 3 on Hubble is already delivering novel and interesting results 
(Berta et al. 2012, Swain et al. 2012). 

Whilst in previous years the transit technique, combined with the photometric precision 
of Hubble and Spitzer, has been an asset for the success of this field, debates have originated 
in the community from the impossibility of repeating observations beyond any reasonable 
doubt. Some of these concerns have been addressed by adopting more robust and objective 
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statistical techniques to remove instrumental systematics (Gregory 2011, Gibson et al. 
2012a, Waldmann 2012, Waldmann et al. 2012a). Furthermore there has been a rapid 
escalation of successful results from the ground over the last few years (Redfield et al. 2008, 
Snellen et al. 2008, 2010, Bean et al. 2010, Swain et al. 2010, Waldmann et al. 2012a, 
Brogi et al. 2012). Ground-based observations have the non-trivial limitation of having 
telluric line contamination interfering with measurements, especially in the infrared where 
most of the key molecules show stronger absorption features. However, observations from 
the ground can be repeated more easily (Waldmann et al. 2012a) and in some cases they 
can cover spectral regions not reachable from space (Swain et al. 2010), or provide higher 
spectral resolution (Snellen et al. 2010, Brogi et al. 2012). 

The hot Jupiter HD- 189733b has been the most observed planet to date due to the 
brightness of its mother star and its favourable atmosphere. Observations, taken with 
multiple instruments and by different teams (e.g., Knutson et al. 2007, Tinetti et al. 2007, 
Charbonneau et al. 2008, Grillmair et al. 2008, Swain et al. 2008, 2009), have been 
interpreted with the presence of H 2 0, CH 4 , C0 2 and CO in its atmosphere (Barman 2008, 
Swain et al. 2009, Madhusudhan & Seager, 2009, Line et al. 2010, Tinetti et al. 2010a, 
Lee et al. 2012). Gibson et al. 2012 and Pont et al. 2012 suggest instead that hazes 
extend from the UV down to the MIR at the terminator. The most recent detection of 
non-thermal emission from its dayside, related to the 3.3 //m methane feature (Swain et al. 
2010, Waldmann et al. 2012a), has increased the interest towards this object even further 
(Mandell et al. 2011). 

Here we present new observations of HD-189733b in the J-, H- and K-bands with the 
NASA Infrared Telescope Facility (IRTF)/SpeX instrument. While only two nights out 
of three could be used for our analysis, they were good enough to allow the extraction of 
low resolution spectra in those bands. Given that the J- , H- and K-bands were observed 
simultaneously, these measurements have a clear advantage for atmospheric interpretation 
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over photometric data recorded at different times (Sing et al. 2009, Gibson et al., 2012b). 
The activity of the star (Knutson et al. 2007, 2009, Agol et al, 2010), in fact, may prevent 
the level of accuracy needed to detect molecular features in the planetary atmosphere 
(Ballerini et al. 2012). Furthermore, most of the current instruments from space and 
ground do not have an absolute calibration at the 10" 4 level. 

To analyse the data, we have used the techniques described in Swain et al. (2010) 
and Waldmann et al., (2012a) tailored for transmission spectroscopy. These techniques 
have been validated on observations with a reference star in the field over the 0.28-4.2yum 
spectral range (L-band), where the telluric contamination is high. Waldmann et al. (2012a) 
found that the spectrum of the observed reference star was flat, indicating an adequate 
background subtraction. 



2. Data reduction and spectral extraction 
2.1. Observations 

We describe here the observations obtained on 12 June 2008 and 22 July 2008, using 
the SpeX spectrograph at the NASA Infrared Telescope Facility. The instrument was 
configured in the short- wavelength cross-dispersed mode (SXD: 0.8-2.5 /iin), in which five 
orders of the spectrograph covered the spectral range from 0.8 to 2.5/iin simultaneously, 
with a small gap in wavelength coverage due to a strong telluric absorption feature between 
1.81 and 1.86/iin (Rayner et al., 2003). We used a slit of l''6xl5'', an average integration 
time of 6 sec and the ABBA nodding sequence. Flat field and argon lamp calibrations were 
obtained before and after the transit. The observations were scheduled to start two hours 
before the primary transit and to stop two hours after the planetary egress. 



-6 - 



Comparing the two nights, 22 July 2008 presents a higher scatter and a gap in the 
lightcurves between -0.032 and -0.028 of the normalized orbital phase, due to a switch-off of 
the instrument. Moreover its signal-to-noise ratio (SNR) is lower than the SNR obtained 
over the night of June 12th (Fig. [1]). An additional transit event has been observed with 
the same settings on 23 June 2010, but was discarded due to the unrecoverable change of 
throughput during the planet ingress. 



2.2. Data Analysis 

We used the SpeX data reduction package, Spextool (Cushing et al. 2004) and we 
applied the standard calibration procedure, which includes background subtraction, flat 
fielding and wavelength calibration. This procedure yielded sets of 382 and 450 individual 
stellar spectra for the first and the second night respectively. We removed the outliers by 
sigma clipping the time-series T\ at 4a and we replaced the bad pixels with the mean of 
the surrounding pixels. We applied a first-order airmass correction, using the cosine-based 
approximation provided by David Tholen in the IDL library. We then obtained the modified 
time-series T' x = T x x ce -«^rm^. 

To extract the planetary signal from the raw data, we had to use de-noising techniques 
followed by amplification of the signal. We thus applied the Model- Correlate- Fit (MCF) 
method (Swain et al. 2010, Waldmann et al. 2012a) to the data, which is useful when 
the wavelength correlated noise is stronger than the channel-to-channel differential signal. 
The MCF approach allows the extraction of the exoplanetary signal by correlating multiple 
time-series in the Fourier domain. Since here we treat primary transit observations, we 
removed the average transit depth in the J-, H- and K- bands by dividing by the mean flux 
of the analysed band. The final spectrum is thus a differential transmission spectrum. 
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A critical step in our analysis is the cleaning of all the non-random temporal changes 
in the spectro-photometric time-series that are not strictly associated with the exoplanetary 
transit (i.e., variations of pixel sensitivity, bias off sets). We therefore normalised each 
spectrum Fi (i corresponds to each spectral frame of the dataset) to remove all the 
systematic errors correlated in wavelength: 

F*(\) = (1) 

where F* is the normalized spectrum, Fi is the flux of each spectrum i as a function of 
wavelength A and Fb is the mean flux of the analysed band, which ranges from Ai = 0.94 
jum to A2 = 1.39 /im for the J-band, from Ai = 1.41/mi to A2 = 1.81//m for the H-band and 
from Ai = 1.94/mi to A2 = 2.42^um for the K-band : 

_ [ A A2 F,(A)dA 

F B = J \ ! (2) 
A2 — Ai 

It is important to stress that the normalisation of the broad-band flux does not erase the 
wavelength-dependent spectral modulations. 

Among the data sets analysed, the one recorded on 12 June 2008 presents fewer 
systematics and is rather uniform in terms of scatter and flux modulations. In contrast, 22 
July 2008 data exhibit a higher scatter, especially during the pre-ingress, most likely due 
to the poorer atmospheric conditions (i.e. probably the presence of a thin layer of cirrus). 
We noticed that for highly scattered time-series, the removal of the mean broad-band flux 
caused an offset between the A and B frames obtained with the nodding technique. These 
problematic time-series, either at the edges of the bands analysed or associated with low 
SNR (Fig. [1]), were discarded to prevent the addition of systematics. One needs to be 
cautious about this process though: while a proper correction must be applied to eliminate 
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the systematics, too few spectral channels do not allow us to extract a meaningful signal. 
The number of data points and the uncertainties in the final spectra (Fig. [TUl fTT]) are 
directly related to this critical step. 

We applied a second normalisation to remove residual fluctuations in the time-series, 
obtaining F**(X). 



F**(\) = 5^ (3) 



a = (Ai + k ■ X), b=(Xi + (K + l)-X) (5) 

where X is the number of spectral channels we combined and Fx is the mean flux of the 
X channels. X should be optimised on the quality of the dataset. We found that usually 
~ 100 to 150 was a good compromise between the SNR requirements and the spectral 
resolution, k represents a set of X channels. The second panel in Fig. [2] shows a lightcurve 
for a single channel after all the corrections. 

After the pre-cleaning process, each lightcurve T" is still too noisy to show the 
transit. We extracted the signal common to all time-series by Fourier-transforming all the 
time-series in the spectral band considered and by stacking together the X-transformed 
spectral channels (Fig. [21 third panel). As a result, the spectral resolution is reduced. 
Finally, we normalised the correlated time-series by taking the geometric mean and the 
inverse Fourier transform: 
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x 

nit) = &-\\[&[n{t)]) 1/x (6) 

A=l 

where & denotes the Fourier transform and J^ -1 its inverse. 

We fitted each of the final lightcurves with two components: a residual baseline curve 
(characterised by a second order polynomial function) fitted to the out-of-transit and a 
transit lightcurve model (Mandel & Agol, 2002) with the transit depth being the only free 
parameter (Fig. [2l bottom panel). We used the limb-darkening coefficients by Claret, 
(2000). The transit depths 5 K of these final lightcurves are the differential transmission 
spectra. We repeated the same procedure for J-, H- and K- bands for both nights. 

To calculate the error bars, we considered both the standard deviation of the model 
substracted residual, (Jscat, and the standard deviation of the transit depth derived from 
changing baseline fits, a fit- The final standard deviation a is obtained by the quadrature 
of the two previous terms: a = \/ ^scat + °\it 

To test the robustness of our results, we used a Markov-chain Monte Carlo method 
(MCMC) to re-fit the lightcurves. We did two MCMC fits, the first one to model the 
residual baseline curve and the second one to fit the transit lightcurves. The posteriors 
of the first fit were used to create a covariance matrix. From the same covariance matrix 
we generated a multivariate prior for the second fit. The retrieved depths and standard 
deviations are consistent with our previous analysis (Fig. [5]). 

Additionally, combining the spectra obtained for the 12 June 2008 night we obtained 
a transmission spectrum of HD-189733b that covers the wavelength range from 0.94 to 2.4 
/im. As the wavelength variation shown by the result is indicative of the presence of spectral 
features, we proceeded to analyse this aspect. We thus compared the final transmission 
spectrum with syntethic models of the planetary atmosphere (see Section HI Fig. [121 n~3j) . 
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Note that this comparison is possible as all the J-,H- and K- band spectra were recorded 
simultaneously. 

3. Models 

The transmission spectra were modeled using a line-by-line radiative transfer code 
as described in Tinetti et al. (2007) and Hollis et al. (submitted). For all models, an 
isothermal atmospheric profiles at T ~ 1500K, probing down to 10 bar, was used. The 
atmosphere was taken to be cloud-free for a start, but Rayleigh scattering due to molecular 
hydrogen (Liou 2002) and H2 - H2 collision-induced absorption (hereafter: CIA - Borysow 
et al. 1997) were included. The CIA continuum (plus Rayleigh scattering at the shorter 
wavelengths) provides an effective absorption floor, on top of which was added absorption 
from H 2 0, CH 4 , CO2 and CO, these being the molecules previously detected on the dayside 
(Swain et al. 2009, Madhusudhan & Seager 2009, Line et al. 2010, Lee et al 2012). The 
mixing ratios (abundances) of these molecules were varied in order to determine the best 
composition, i.e., to provide an adequate fit to the observations. 

For most of the molecules considered, theoretical line lists from the Exomol projectQ 
(Tennyson et al. 2012) were used. For example, for the H2O molecule, absorption 
cross-sections were obtained from Barber et al. (2006), which contributes absorption 
features in the 0.94-2.42 micron range in particular. The CH 4 opacity was modelled using a 
preliminary version of a new variational line list at high temperature (T=900K; Yurchenko 
et al. in prep.). For the CO and CO2 molecules, the HITEMP line list (Rothman et al. 2010) 
was used. However at wavelengths less than A = 1.43 /im these data were not available, 
and data from hitran line list (Rothman et al. 2009) at T = 300K were used. 



1 http: / / www.exomol.com 
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4. 



Results 



Using the method described in the previous section, we analysed the J-,H- and K-band 
data for both nights. By combining sets of 150 channels we extracted the differential depths 
as a function of A (Tabled]). We show in Figs. Eland [7] the lightcurves for the first night, J- 
and H-, K- respectively. Notice that this approach does not provide an accurate absolute 
calibration, but rather a refined relative measurement. Uncertainties much higher than 
the photon-noise and the residual correlation are typical of ground-based observations and 
depend on the quality of the data. The uncertainties in our data are, on average for the first 
night, three times larger than photon noise in the H-band and two times larger in the J- 
and K-band. For the second night the uncertainties for the H- and K- band are three times 
larger than the photon noise. The J-band on the second night showed a lower signal level. 

To test the statistical significance of our results we estimated the auto-correlation 
function (ACF) for our time-series before and after the data reduction (Fig. [8]). Given a 
discrete signal S n , the ACF(r) is the cross-correlation of the signal S n with itself, at lag r. 



The auto-correlation function is a mathematical function that, over consecutive time 
intervals, gives the degree of similarity between a given time-series and its lagged version. 
In a photon-noise limited case the ACF of a lightcurve will show a peak at r = and will 
have zero value over the other lags. In Fig. |H]we present the ACF over 370 lags for one 
raw lightcurve and for the residual to the final lightcurve. The raw data present a strong 
correlation. By contrast, the autocorrelation for final residuals is all within la confidence 
and their amplitudes negligible. This indicates that correlations in time were efficiently 




(7) 



n 
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removed, that our results are normally distributed and white-noise dominated at most 
frequencies. At lower lags, some residual correlation is still persistent but these systematics 
are at frequencies higher than the time scale of the transit event of HD-189733b, hence the 
signal should not be directly affected. 

Due to the really good transmission window within 0.94-1.30/im (Fig. [3]) and the 
good quality of the 12 June 2008 night, we compared the J-band spectrum to synthetic 
models that included H 2 0, CH 4 , CO2 and CO, these being the molecules detected in the 
dayside of the planet (Fig. [T3"|) . We tested one molecule at a time and then combinations of 
these. The best fit for the J band was obtained with water with a mixing ratio of 5-10 -4 at 
temperature T ~ 1500K. We then compared the best model achieved with the whole 0.94 

- 2.42/iin spectral range. Note that our J-, H- and K- spectra are relative measurements, 
hence their absolute level is unconstrained. It is therefore possible to renormalise their 
offset. Thus, keeping our J-band spectrum fixed, we matched the H- and K-band spectra 
to the model, finding an excellent agreement. The best fit was given by the model that 
includes water with a mixing ratio of 5-10 -4 at a temperature that ranges between T~1000 

- 1500 K. We tested as well the possibility of a flat spectrum by fitting a straight line to the 
data but the fit is really poor, as confirmed by the x 2 (see Fig. [12]) . 



4.1. Removal of telluric contamination in the data. 

The entire MCF method is an efficient procedure to remove, through several steps, 
the telluric contamination contained in the raw data. The raw data contain both the 
signal and the noise components. These components are difficult to disentangle in the 
time-domain, but can be described by different sets of frequencies in the Fourier domain. 
This is why, when we combine together ~ X fourier-transformed lightcurves, noise such as 
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the telluric contamination is weakened, while the transit signal is strengthened. We refer to 
Waldmann et al. (2012a) for a detailed explanation. To quantify the residual contamination 
by systematics we applied the procedure described in Sect. 12.21 only to the out-of-transit 
lightcurve, i.e., removing the eclipse signal. We find that the amplitude of the systematic 
noise and the residual telluric component is within the uncertainties of the planetary signal 
over the whole spectrum. The effectiveness of the MCF was tested with the ACF. 
Additionally, we report the weather condition over the two nights. Figure H] shows 
the temperature, relative humidity and pressure readings of the Canada-France-Hawaii 
Telescope (CFHT)g weather station as well as the atmospheric opacity at 225 GHz by 
the Caltech Submillimeter Observatory (CSoJf. No significant correlations between these 
parameters and the expected transit shape was found. 

We show in Fig. [3] the atmospheric trasmission on Mauna Kea, obtained in similar 
conditions to our nights (UKIEr^ web pages). The intensity of the transmission is on 
average 1 = 1, with some absorption features around 1.125/mi, 1.41/xm, and 2.06/im whose 
effect reflects directly on the error bars of the data points evaluated at the same A. Those 
error bars are, in fact, the largest in the J-, H- and K- band respectively. The strongest 
telluric absorptions are reported around 1.355/im and 2.0/zm; hence the decision to flag the 
two measurements at those wavelengths as less reliable. 



2 http: / / mkwc.ifa.hawaii.edu/archive 
3 http: / / ulu.submm.caltech.edu 

4 http: / /www .jach.hawaii.edu/UKIRT/astronomy/utils/atmos-index.html 
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5. Discussion 

While the spectra obtained during the two nights are consistent with each other within 
the error bars, the spectra extracted on the first night are of superior quality. Among 
the hurdles encountered, we can list the small gap in the flux during the orbital phase, 
the instrumental systematics (e.g., variations of the pixel sensitivity), the residual telluric 
contamination, and the presence of high altitude clouds during the second night. All these 
effects may influence our results at different levels but, by applying the MCF method to 
the IRTF data, we have been able to extract statistically significant planetary spectra. 
The ACF test, in fact, proved that our results are normally distributed and white noise 
dominates at most frequencies. Nonethless these results could be largely improved by 
combining together a large number of good nights. As demonstrated by Waldmann et al. 
(2012a), we can potentially amplify our signal by a factor ~ N by combining N multiple 
good observations. In the work presented here one night was outstanding, one was less 
than average, one was unusable. We decided to flag the measurements at 1.35/mi and 
2/mi (marked in Figs. [121 H3] in yellow), as less reliable. At these wavelengths, the telluric 
contamination is very high and it might affect the measurements (Fig. [3]). On top of that 
the data point at 1.35yum includes fewer spectral channels than our standard procedure, 
since it lies at the end of the band. 

5.1. Comparison with previous observations 

We compared our J-band results with the Hubble /ACS spectrum (Pont et al. 2008) as 
there is a marginal overlap with the two passband points in the 0.95 and 1.05 /im range. 
We find that the results are consistent at the l-er level (Fig. [9]). The comparison with the 
two photometric points recorded with HST/WFC3 (Gibson et al. 2012b) is less obvious, as 
these data points are recorded at different epochs and the star is quite active. Nevertheless 
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they are within the errorbars. 

The H- and K-band spectra have been compared to the transmission spectrum of 
HD-189733b observed by Swain et al. 2008 (hereafter SW08), as it probes the same spectral 
range. If we rebin our measurements to the SW08 spectral resolution we find that our IRTF 
differential transmission spectra are consistent with SW08 at the 2-a level in the H-band 
(Fig. [JO]) and at the l-a level in the K-band (Fig. [TT]) . 

All data are plotted with uncertainties of ± la and centred on the median wavelength of 
each set k of channels. 

Two re-analyses of SW08 were made by Gibson et al. (2012a) and Waldmann et al. (2012b). 
Our data are consistent with both reanalyses at the l-a level. Notice that our errorbars are 
smaller that the ones estimated by Gibson et al. (2012a) but slightly larger than the ones 
by Waldmann et al. (2012b). A detailed comparison between Gibson et al. (2012a) and 
SW08 has been discussed in Waldmann et al. (2012b); we refer to that paper for further 
considerations on this issue. 

Given that the IRTF ground-based measurements and the Hubble space-based 
measurements have completely different sources of systematic errors, the consistency among 
ground and space observations is quite remarkable. 

5.2. Data interpretation 

As mentioned in Section H] we used the 12 June 2008 J-band spectrum and simulated 
models to calibrate the H- and K- spectra recorded during the same night. The models 
generated included the contributions of H2O, CH4, CO2 and CO. The best fit to the data 
was given by the water opacity. The retrieved abundance of water is somewhere in the range 
1 - 5-10 -4 if in the models the floor is due to collision-induced absorption and Rayleigh 
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scattering. If the floor is due to a cloud deck or haze opacity, then this abundance could be 
underestimated. Our data are not sensitive enough to establish the additional contributions 
of molecules such as CH 4 , CO2, CO detected by other teams in the day side of the planet, 
or the terminator. The detection of water vapor is, though, a sound conclusion. To date, 
water vapor absorption in those bands has been reported on two other hot-jupiters, XO-lb 
and XO-2b (Tinetti et al. 2010b, Crouzet et al. 2012) and one super-earth (Berta et 
al. 2012). Hazes have been suggested by Pont et al. (2008) to explain the absence of 
molecular features in the visible transmission spectrum of HD-189733b, yet at the same 
time Redfield et al. (2008) detected the presence of sodium in the same spectral range. 
Gibson et al. (2012b) and Pont et al. (2012) extend the opacity due to hazes down to the 
mid-infrared, proposing that the transmission spectrum is flat from the UV down to 8/rni. 
Our observations do not support this hypothesis in the near infrared, as highlighted by 
Figure [T2l and the % 2 test. Additional work remains to be done to determine the precise 
molecular abundances of water and other species. To this end, observations with broad, 
instantaneous spectral coverage will be especially valuable, perhaps critical, to avoid any 
significant problem associated with stellar variability. 



6. Conclusions 

We have presented here the first ground-based spectroscopic observations of the 
primary transit of the hot- Jupiter HD-189733b, recorded with the NASA IRTF/SpeX 
instrument. We have pre-cleaned our data and applied the Model Correlation Fit technique, 
finding that the J-, H- and K- band spectra are consistent with the collection of datasets 
recorded from space with Hubble. The auto-correlation test demonstrated that correlations 
are efficiently removed, that our residuals are normally distributed and that white-noise 
dominates at most frequencies. By comparing the J-, H- and K- 12 June 2008 spectra to 
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synthetic models, we found that water vapour with a mixing ratio computed between 10~ 4 
and 5-1CT 4 explains the spectral modulations between 0.94 - 2.4/im. Our results are not 
sensitive enough to give further constraints on other molecules like methane, carbon dioxide 
or carbon monoxide detected in other datasets; however by combining together multiple 
good nights the accuracy needed could be easily achieved. 

The work presented in this paper shows that low-resolution exoplanet spectroscopy is 
indeed feasible with medium-sized telescopes from the ground. While the telluric absorption 
is a non-neglegible hurdle to sound exoplanetary atmospheres, the potential to repeat the 
observations with relative ease makes the ground an appealing and complementary option 
to space. 
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Fig. 1. — Signal to noise ratio measured over the totality of the eclipse event for both nights 
of observation. 
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Fig. 2. — From the top: (1) Raw flux for a single pixel-based channel in the J-band, (2) Light 
curve of a single channel after airmass correction and normalization, (3) A sample of three 
sets of 150 channels correlated together in the Fourier domain, (4) Detrended and binned 
light curve of one of the J-band set with an eclipse model overplotted. 
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Fig. 3. — Atmospheric transmission on Mauna Kea over the three analysed bands. These 
data, produced using the program IRTRANS4, were obtained from the UKIRT web pages, 
(airmass = 1.0, H 2 = 1.2 mm, R = 3000) 




Fig. 4.— Weather over 12 June 2008 night (red) and 22 July 2008 (blue), from the CFHT 
Weather station. Top to bottom: temperature (°C), relative humidity (%), pressure (mb) 
and the optical depth, tau (225 GHz) from CSO. The dashed vertical lines mark the transit 
event duration. 
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Fig. 5. — Histograms of the depths retrieved with the MCMC {12 June 2008). From the top: 
lightcurve at 0.958/mi (J-band), lightcurve at 1.459/iin (H-band) and lightcurve at 2.08/im 
(K-band). Note that the standard deviations measured with the MCMC are smaller than 
the ones computed with the MCF technique. 
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Fig. 6. — The 12 June 2008 detrended lightcurves for the J-band showing in red the fitted 
differential lightcurve model. The uncertainties of the individual measurements are ± la 
deviation. For the clarity of the reader we shifted the lightcurves. The wavelengths are 
increasing from bottom to top and left to right. 
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Fig. 7. — The 12 June 2008 detrended lightcurves for the H-band (left) and the K-band 
(right) showing in red the fitted differential lightcurve model. The uncertainties of the 
individual measurements are ± la deviation. For both bands all the lightcurves are shifted 
for clarity and the wavelength is increasing from the bottom to the top. 
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Fig. 8. — Showing the auto-correlation function of one raw lightcurve (black diamonds) and 
one final lightcurve obtained after our data analysis (red points) for the H-band of the 12 
June 2008 night. The 3a, 2a, la confidence limits that the data are normally distributed 
are plotted in blue (solid line), green (dashed line) and red (dash-dot line). All the lags are 
within la limit showing that the correlations were efficiently removed, leaving a white noise 
dominated signal at most frequencies. 
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Fig. 9. — 12 June 2008 J-band differential transmission spectrum of HD-189733b with 
uncertainties of ± la. For comparison the two ACS bandpass points (Pont et al. 2008) are 
shown in red. 
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Fig. 10. — H-band differential transmission spectrum of the hot Jupiter HD-189733b (black 
points), with uncertainties of ± la. For comparison the SW08 spectrum is shown in red. 
(Left: 12 June 2008, Right: 22 July 2008) 
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Fig. 11. — K-band differential transmission spectrum of the hot Jupiter HD-189733b (black 
points), with uncertainties of ± la. For comparison the SW08 spectrum is shown in red. 
Left: 12 June 2008, Right.22 July 2008. 
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Fig. 12. — 0.94-1.4 /im IRTF transmission spectrum compared to a simulated spectrum for 
water with a mixing ratio of 5 • 10 -4 , assuming an isothermal atmosphere at T ~ 1500K. 
The x 2 value for this fit is also given, where the data values known to suffer from telluric 
contamination (marked 'Telluric' in the figure legend) were excluded for this calculation. For 
reference, a straight line with a value equal to the mean of all of the data points (y = 0.0243) 
is also shown, with the associated x 2 value. 



-27- 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (microns) 



H20 (5e-4) 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (microns) 



H20 (5e-4) 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (microns) 



Fig. 13. — 0.94-1.4 /iin IRTF transmission spectrum compared to different families of simu- 
lated spectra. Each family includes water with mixing ratio 5 ■ 10~ 4 plus varying abundances 
of methane (top), carbon dioxide (center) and carbon monoxide (bottom), assuming an 
isothermal atmosphere at T ~ 1500K. Water alone can explain well the features of our IRTF 
transmission spectrum; the data are not sensitive enough to give constraints on the presence 
of other molecules. 
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2.392 
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Table 1: Differential depth values and relative errors for 12 June 2008 night (left) and for 
22 July 2008 night (right). All the wavelengths are in microns. The telluric points are not 
listed. 
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